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Axions currently provide the most compelling solution to the strong CP problem. These particles 
may be copiously produced in the early universe, including via thermal processes. Therefore, relic 
axions constitute a hot dark matter component and their masses are strongly degenerate with those 
of the three active neutrinos, as they leave identical signatures in the different cosmological observ¬ 
ables. In addition, thermal axions, while still relativistic states, also contribute to the relativistic 
degrees of freedom, parameterised via N e g. We present the cosmological bounds on the relic axion 
and neutrino masses, exploiting the full Planck mission data, which include polarization measure¬ 
ments. In the mixed hot dark matter scenario explored here, we find the tightest and more robust 
constraint to date on the sum of the three active neutrino masses, < 0.136 eV at 95% CL, 

obtained in the well-known linear perturbation regime. The Planck Sunyaev-Zeldovich cluster num¬ 
ber count data further tightens this bound, providing a 95% CL upper limit of 'Y)rn l , < 0.126 eV 
in this very same mixed hot dark matter model, a value which is very close to the expectations in 
the inverted hierarchical neutrino mass scenario. Using this same combination of data sets we find 
the most stringent bound to date on the thermal axion mass, m a < 0.529 eV at 95% CL. 

PACS numbers: 98.80.-k 95.85.Sz, 98.70.Vc, 98.80.Cq 


Introduction — The axion field arises as a solution to 
solve the strong CP problem in Quantum Chromodynam- 
ics m- The axion is the Pseudo-Nambu-Goldstone as¬ 
sociated with a new global U(1)pq (Peccei-Quinn) sym¬ 
metry that is spontaneously broken at an energy scale 
f a . In the early universe, axions can be produced via 
thermal or non-thermal processes. While in the former 
case, the axion contributes an extra hot thermal relic (to¬ 
gether with three active sterile neutrinos), in the latter 
case, the axion could be the cold dark matter component. 
In the following, we shall focus on the thermal axion sce¬ 
nario. In order to compute the present thermal axion 
relic density, the most relevant process is the axion-pion 
interaction, 7r + 7r — > n + a. The characteristic parameter 
for the thermal axion is / a , the axion coupling constant, 
that can be related to the axion mass by 


/ttTOtt VR ,, 10 7 GeV 

m - = — rrs =0 - 6eV -^ 


di 


where the up-to-down quark mass ratio is taken as R = 
0.553 ± 0.043, and = 93 MeV is the pion decay con¬ 
stant. 

Thermal axions, while still relativistic, will increase 
the amount of radiation in the universe, contributing to 
the effective number of relativistic degrees of freedom 
N e s. In the standard cosmological ACDM model with 
three active neutrino species, we expect N e g = 3.046 [4], 
where the 0.046 takes into account corrections for the 


non-instantaneous neutrino decoupling from the primor¬ 
dial plasma. An extra AN e g = N e g — 3.046 modifies 
the damping tail of the Cosmic Microwave Background 
(CMB) temperature angular power spectrum, changing 
two important scales at recombination, the sound hori¬ 
zon and the Silk damping scales, as well as the primor¬ 
dial abundances of the light elements predicted by Big 
Bang Nucleosynthesis. When thermal axions become 
non-relativistic particles, they will affect the different cos¬ 
mological observables in an analogous way to that of mas¬ 
sive neutrinos, i.e. by increasing the amount of the (hot) 
dark matter density in our universe. Axions will sup¬ 
press structure formation at scales smaller than the free- 
streaming scale, favouring clustering only at large scales. 
Thermal axions will also leave an imprint on the CMB 
temperature anisotropies, via the early integrated Sachs- 
Wolfe effect. Therefore, a large degeneracy between the 
axion mass and the total neutrino mass is expected [5]. 
Several papers in the literature have provided cosmolog¬ 
ical constraints on the thermal axion mass in different 
cosmological scenarios, see e.g. Refs. tSHIB. 

In light of the recent Planck 2015 temperature and po¬ 
larization data [12j . it is timely to compute the changes 
in the existing bounds on the thermal axion mass, in¬ 
cluding the case in which massive neutrinos are present. 
Our results are obtained using the Monte Carlo Markov 
Chains (MCMC) package CosmoMC [13], with CAMB (Code 
for Anisotropies in the Microwave Background) [IT] as 
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Parameter 

Prior 

fib h A 

[0.005,0.1] 

^edm^ 

[0.001,0.99] 

e s 

[0.5,10] 

T 

[0.01,0.8] 

m a (eV) 

[0.1,3] 

E (eV) 

[0.06,3] 

n s 

[0.9,1.1] 

log[10 10 A s ] 

[2.7,4] 


TABLE I: Priors for the parameters used in the MCMC anal¬ 
yses. 

a solver of the Boltzmann equations. In the mixed hot 
dark matter scenario, in which both axion and neutrino 
masses are allowed to freely vary, we find the tightest and 
most robust constraint to date on the sum of the three 
active neutrino masses, E m i' < 0.136 eV at 95% CL, as 
this only relies on the (very well-known) linear perturba¬ 
tion regime. 

Thermal axion cosmological model The scenario 
we analyze here is the ACDM model, with both axions 
and neutrinos as extra hot thermal relics. We describe 
this scenario by the following set of parameters: 

{w 6 ,w c ,0 s ,r, m a , E «V,n s ,log[10 10 A s ]} , (2) 


of Ref. [T5]. We also include the new Planck 2015 lens- 
ing likelihood EH, constructed from measurements of the 
power spectrum of the lensing potential, referring to it 
as lensing. Concerning Planck catalogs, we make use of 
the Sunyaev-Zeldovich second cluster catalog [TJ. TTi| (de¬ 
noted as SZ in what follows), which consists of 439 clus¬ 
ters with their corresponding redshifts and with a signal- 
to-noise ratio q > 6. We also consider additional datasets 
to the Planck satellite measurements, as a gaussian prior 
on the Hubble constant Hq = 73.8 ± 2.4 km/s/Mpc, ac¬ 
cording with the measurements of the Hubble Space Tele¬ 
scope, [2D] - We refer to this data set as HST. We also 
include measurements of the large scale structure of the 
universe in their geometrical form, i.e. in the form of 
Baryon Acoustic Oscillations (BAO). In particular, we 
use the 6dFGS, SDSS-MGS and BOSS DR11 measure¬ 
ments of Dy/ r 2 ^ [2TK23] . referring to the combination of 
all of them as BAO. We shall also consider large scale 
structure measurements in their full matter power spec¬ 
trum form, as provided by WiggleZ survey [23], and de¬ 
noted as MPK. Tomographic weak lensing surveys pro¬ 
vide a powerful tool to constrain the mass distribution 
in the universe, and therefore we shall also exploit in our 
analyses the constraint on the relationship between ag 
and Q m of cr8(f2 m /0.27) 0 - 46 = 0.774 ± 0.040 provided by 
the Canada-France-Hawaii Telescope [25], CFHTLenS. 
This last measurement is referred to as WL. 


where uib = flbh 2 is the baryon matter energy density, 
uj c = f l c h 2 the cold dark matter energy density, Q s is the 
ratio between the sound horizon and the angular diam¬ 
eter distance at decoupling, r is the reionization optical 
depth, m a is the axion mass in eV and E m„ the sum of 
three active neutrino masses in eV. We consider also the 
inflationary parameters, the scalar spectral index n s and 
the amplitude of the primordial spectrum A s . We use flat 
priors for all the parameters, as listed in Tab. |T] Notice 
that the standard extra radiation density will change, as 
the presence of a thermal axion will increase the value of 
the effective number of relativistic degrees of freedom in 
the following way: 


A-ZVeff 


4 /3 na\ 4/3 
7 \ 2 n v ) 


(3) 


where n a is the axion number density and n„ is the 
present neutrino plus antineutrino number density per 
flavor. The current axion number density is a function 
of the axion decoupling temperature To, that is a func¬ 
tion of the axion mass m a . For the details related to the 
calculation of the axion decoupling temperature, we refer 
the reader to Ref. m- 

Datasets Our baseline data set consists of the recent 
Planck 2015 satellite CMB temperature and polarization 
measurements mmm- We consider a combination of 
the likelihood at 30 < £ < 2500 using TT, TE and EE 
power spectra and the Planck low-£ multipole likelihood 
in the range 2 < i < 29. We refer to this combination 
as Planck TT,TE,EE+lowP , following the nomenclature 


Results Table [TT] summarises the results from our 
MCMC analyses in the mixed hot dark matter scenario 
revisited here. Notice that Planck temperature and po¬ 
larization measurements ( TT,TE,EE+lowP ) set 95% CL 
upper bounds of E m v < 0.441 eV and m a < 2.09 eV 
respectively. The bounds on the thermal axion mass are 
similar to those obtained in the case in which only ax¬ 
ion masses are considered, albeit for that case the value 
of the (is parameter is always higher than the one shown 
here, as only one hot relic suppresses the small-scale clus¬ 
tering. Nevertheless the deviation of as is not significant 
(about a half-sigma away from the value illustrated in 
Tab. 0- Furthermore, neutrino oscillation experiments 
have provided compelling evidence for the existence of 
neutrino masses and therefore neutrinos must be added 
as massive particles. The addition of CMB lensing mea¬ 
surements from the Planck satellite weakens the neutrino 
mass bounds, as discussed in m ■■ the lensing reconstruc¬ 
tion data prefers lensing amplitudes lower than the stan¬ 
dard prediction, and this favours higher neutrino masses, 
as the presence of these will smooth the lensing power 
spectrum. Summarizing, when Planck CMB lensing con¬ 
straints are considered, the neutrino mass bounds are 
pulled away from zero, and we obtain E m„ < 0.538 eV 
and m a < 1.67 eV at 95% CL. The addition of weak 
lensing constraints on the relationship between the mat¬ 
ter clustering amplitude 0 $ and the matter mass-energy 
density 0 TO to Planck TT,TE,EE+lowP measurements 
only mildly tightens both the thermal neutrino and ax¬ 
ion masses. The largest impact on both E m v and m a 
bounds comes from the large-scale structure information 
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TT,TE,EE+lowP 

TT,TE,EE+lowP 

+lensing 

TT,TE,EE+lowP 

+WL 

TT,TE,EE+lowP 

+MPK 

TT,TE,EE+lowP 

+BAO 

TT,TE,EE+lowP 

+HST 

TT,TE,EE+lowP 
+BAO +HST 

TT,TE,EE+lowP 
+BAO +HST +SZ 

n c ft 2 

n 1 oqk+0-0034 
u.x^oo_q Q036 

0 1 9QK+0.0034 

u.xz,oo_o 0034 

n 1 ook+ 0.0032 
o Q032 

fl 1OQ7+0-0034 

v - lzo ' -0.0031 

O 1 990 + 0 .0023 

u.izzo_o 0023 

O 1 99q+0.0032 

u.±zz«3_o 0032 

0 1220" 1 " 0 ' 0024 
u.±z,z,u_o .0023 

0 1216’*’ 0-0023 
u.xzxu_o .0023 

m a [eV] 

< 2.09 

< 1.67 

< 1.87 

< 0.835 

< 0.763 

< 1.21 

< 0.709 

< 0.529 

£% [eV] 

< 0.441 

< 0.538 

< 0.360 

< 0.291 

< 0.159 

< 0.182 

< 0.136 

< 0.126 

O' 8 

0 779 + 0 - 083 
u -' ‘ y —0.094 

0 767+ 0 - 065 
u - -0.072 

0-7891“;“™ 

D R14+0-049 

U.814—o 056 

0.827l»;°42 

0.8201“;“;) 

0-829l“;°39 

fl QqK+0-033 

U.O-3D —0 035 

Q m 

0 342+ 0 - 054 
u.o^_q 048 

0 344+°-°55 
0.048 

0-328t°;°^ 

0 326"* -0 ' 033 
u.ozu _ 0 029 

fl “31 9+0-016 
u.oiz,_ 0 .014 

0.3l5l“;“ 3 2 ) 

0.3091“;““ 

0 306’*’ 0 014 
U.- 3 UO -0 013 

Iog[10 10 A s ] 

o lcn +0.067 
-3.l-3l_o.o70 

3 109"*" 0 064 
o.iuy_o 062 

Q n 7+0.071 
°- 11 1 -0.068 

o 19q +0.066 

O-l^l_0.071 

0 1 9«+0 066 

-3.1ZO_Q.070 

q 1 9Q+0.066 

-3.izy_Q.068 

q 1 ^Q+0.065 
-3.1ZO_0 069 

q -■ q O +0.063 

-3.l-5Z_o.064 

n s 

0 972"*" 0-011 
u.y iz_o.oi 2 

0 972’*" 0-010 

u.y iz_o.oii 

0 Q74 '*' 0-011 
u-y <4_o.oi2 

0.97278t°;™9 

0 Q7‘ q 4+ 0 - 0093 

U.y <04_o 0089 

0 Q76+0-0 10 

u-y <o_o.oio 

0 q7fjQ+0-0095 

u.y <09_o .0091 

0.9768l“;™89 


TABLE II: 95% CL constraints on the parameters of the mixed hot dark matter scenario explored here (the ACDM+m a +Sm„ 
model) for the different combinations of cosmological data sets. 


as well as from the prior on Hq from the HST experiment. 
Notice that the bounds are significantly tighter when one 
of the former constraints are considered in the analyses. 
Concerning the Hq prior, the 95% CL upper bounds on 
the thermal relic masses become YZ 171 " < 0.182 eV and 
m a < 1.21 eV. The reason for this large improvement is 
due to the large degeneracy between and H 0 [2S| • 

When Y2 m v increases there is a shift in the distance to 
last scattering. This shift can be easily compensated by 
lowering Hq, resulting in a strong degeneracy between 
these two parameters, which can be broken via an in¬ 
dependent measurement of Ho. However, the tightest 
axion and neutrino mass constraints arise when large- 
scale structure data is exploited in its geometrical form, 
via the BAO signature. Indeed, it was shown in Ref. m 
that, when constraining Y2 in minimal schemes as the 
one explored here (i.e. a ACDM model), the informa¬ 
tion contained in the broadband shape of the halo power 
spectrum was superseded by geometric information de¬ 
rived from the BAO signature. We find here a similar 
effect, although the BAO measurements that we exploit 
correspond to several redshifts and surveys, while the 
full-shape data come from only one survey, the WiggleZ 
survey. Using the full matter power spectrum measure¬ 
ments from the former experiment, we obtain 95% CL 
upper bounds of YZ m v < 0.291 eV and m a < 0.835 eV. 
The 95% CL upper bound of YZ m v < 0.159 eV for the 
Planck TT,TE,EE+lowP and BAO combination is very 
close to the one quoted by the Planck collaboration for 
the same data sets, YZ m v < 0-17 eV [15] . However, 
our constraint is tighter, as we are also considering here 
axions as additional thermal relies, and there exists a 
strong degeneracy among these YZ m u and m a . Figure [l] 
illustrates such a degeneracy. We depict, in the £m„, 
to 0 ) plane, the 68% and 95% CL contours arising from 
the analyses of Planck TT,TE,EE+lowP data plus ad¬ 
ditional measurements, as the Planck lensing signal and 
other data sets ( WL. BA O, HST and SZ cluster number 
counts). Notice that the constraints are greatly improved 
for the former two cases, leading to very tight constraints 


on the masses of these two thermal relics. 

The addition of the BAO and the HST data sets leads 
to the strongest constraint on the neutrino mass to date 
in the linear perturbation regime, YZ m » < 0.136 eV at 
95% CL. The corresponding bound on the axion mass 
is m a < 0.709 eV. The authors of [28] have recently re¬ 
ported, using the one-dimensional Lyman-a forest power 
spectrum of the BOSS experiment, a 95% CL upper 
bound of Y2 m v < 0.12 eV in the case in which only 
massive neutrinos are present. Notice however that this 
constraint strongly relies on hydrodynamical simulations, 
while our bounds are derived in the very well-known lin¬ 
ear perturbation regime. Furthermore, the addition of 
the Planck SZ cluster number counts data provide a 
competitive 95% CL upper limit of YZ 7 * 1 " < 0.126 eV 
in the mixed axion-neutrino hot dark matter scenario 
(the corresponding bound on the thermal axion mass is 
m 0 < 0.529 eV). This limit is very close to the expecta¬ 
tions for Y2 fn v in the inverted hierarchical neutrino mass 
scenario, highlighting the fact that improved cluster mass 
calibrations could help enormously in disentangling the 
neutrino mass spectrum. 

Conclusions The polarization measurements from 
the Planck 2015 data release offer a unique opportunity 
for testing the dark matter paradigm. These recent re¬ 
sults point to a standard ACDM as the preferred model 
for the universe we observe today. Nevertheless, a small 
hot dark matter component can still be present. We 
have explored the most general scenario, i.e. a mixed 
hot dark matter model with two thermal relics, neutri¬ 
nos and axions, which would account for the small con¬ 
tribution from the hot dark matter sector to the total 
mass-energy density of the universe. Using Planck tem¬ 
perature and polarization data, and making use of the 
Planck Sunyaev-Zeldovich cluster catalog as well as in¬ 
dependent, low redshift probes, including measurements 
of the Baryon Acoustic peak in galaxy clustering and of 
the Hubble constant, we derive the tightest bounds to 
date on the thermal relic masses. The 95% upper lim¬ 
its extracted from the numerical analyses carried out in 
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FIG. 1: 68% and 95% CL allowed regions in the (J3 m m a ) 
plane, both in eV, for some of the cosmological data combi¬ 
nations explored in this analysis. 


this study are m a < 0.529 eV and 53 m v < 0.126 eV for 
the axion and total neutrino mass, respectively. These 
results strongly motivate the need for improved cluster 
mass calibrations. They also clearly illustrate the power 
of combining low and high redshift probes when corner¬ 
ing the dark matter thermal properties. 
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